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ABSTRACT 

Some recent investigations of spectropolarimetric observations of the Zeeman effect in the Fe I lines 
at 630 nm carried out with the Hinode solar space telescope have concluded that the strength of the 
magnetic field vector in the internetwork regions of the quiet Sun is in the hG regime and that its 
inclination is predominantly horizontal. We critically reconsider the analysis of such observations and 
carry out a complete Bayesian analysis with the aim of extracting as much information as possible 
from them, including error bars. We apply the recently developed Bayes-ME code that carries out a 
complete Bayesian inference for Milnc-Eddington atmospheres. The sampling of the posterior distri- 
bution function is obtained with a Markov Chain Monte Carlo scheme and the marginal distributions 
are analyzed in detail. The Kullback-Leibler divergence is used to study the extent to which the 
observations introduce new information in the inference process resulting in sufficiently constrained 
parameters. Our analysis clearly shows that only upper limits to the magnetic field strength can be 
inferred with fields in the kG regime completely discarded. Furthermore, the noise level present in the 
analyzed Hinode observations induces a substantial loss of information for constraining the azimuth 
of the magnetic field. Concerning the inclination of the field, we demonstrate that some information 
is available to constrain it for those pixels with the largest polarimetric signal. The results also point 
out that the field in pixels with small polarimetric signals can be nicely reproduced in terms of a 
quasi-isotropic distribution. 

Subject headings: magnetic fields — Sun: atmosphere - Sun: magnetic fields — line: profiles - 
polarization 
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1. INTRODUCTION 

Particularly during the last decade, we have witnessed 
an increasing interest in the investigation of the mag- 
netism of the quiet solar photosphere. There are at least 
two reasons for this. On the one hand, improvements in 
the sensitivity of spectro-polarimeters have allowed us to 
carry out spectropolarimetric observations of great qual- 
ity in regions of very low magnetic flux. On the other 
hand, since the quiet Sun covers a fraction larger than 
90% of the solar surface, it is important to investigate 
whether the magnetism of the quiet photosphere plays 
a significant role on the heating of the outer regions. In 
case the magnetic energy stored in these apparently non- 
magnetic zones turns out to be substantial, as claimed by 
iTruiillo Bueno et alj (|2004f ) , its impact on the overall en- 
ergy balance of the solar atmosphere should be accounted 
for. 

The Zeeman polarization signals observed in the 
quiet Sun are very weak and, even with the best 
spectropolarimeters, it turns out difficult to detect 
them and to interpret reliably the observations. One 
of the main causes for such weak polarization signals 
is that the filling factor of the magnetic field in the 
resolution elements of ground base d telescopes is yer> 
small , of the order of 1% (e.g.. IStenflol [iMQ pr 



19951 : iDommguez C erdena ct al. 2003; K homenko et al 
20031 : IMartmez Gonzalez et alJ l2006al ). The 



reason 



for this small filling factor may lie in the fact that 
the field is organized at very small scales and Stokes 
V Zeeman signals cancel effectively when averaging 



fields with several inclinations in the resolution element 
(linear polarization also cancels when the azimuth is 
also random). Understanding the probability distri- 
bution of field strengths h as p r oduced a long debate 
that is n ot still resolved dLinl Il995t iLin fe Rimmel ' 
19991: IS ocas- Navarro fe Sanchez Almeidal 1200 
Dommguez Cerdena et ail 120031: iKhomenko et al 



tie 
ele 
32; 



2003; ISocas- Navarro et alj l2004t iLites fe Socas-Navarr. 



Martinez Gonzalez et al 



2006bllal: 



Dommguez Cerdena et al.l 120061: lOrozco Suarez et al.l 
2007d: IMartmez Gonzalez et al.1 12008b: 

Ramirez Velez et al.l [2008). The reason is that di- 
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agnostic tools based on the Zeeman effect suffer from 
cancellation effects and it is not clear to what extent 
the probability distribution inferred from the residual 
Zeeman signals can be compared to the probability 
distribution of the field strength inside the pixel. 
Therefore, it is important to develop and apply alter- 
native diagnostic tools that do not suffer from such 
cancellations. In th i s respect, w e note the efforts of 
iTrujillo Bueno et al.l (|2004 120061 ) applying techniques 
based on the Hanle effect to infer properties about 
the magnetic field in the quiet Sun. They concluded 
that there is a substantial amount of "hidden" mag- 
netic energy in the internetwork regions of the quiet 
photosphere. Most of these "hidden" fields remain 
undetected to the Zeeman effect because they are tan- 
gled at scales below the resolution of present telescopes 
or because they are inherently weak. Others have 
applied Zeeman-based techniques in lines with strong 
hyperfine perturbations, which mak e them sensitive to 
the strength of the magnetic field (|L6pez Ariste et all 



Fig. 1. — The left panel shows the observed intensity (normalized 
Fe I doublet. The background image on the right panel presents the 
V 2 ) 1 / 2 d\/ (I) . The red and blue contours delineate regions of f(Q 2 
The white contours indicate regions where I C /(I C ) = 1.02. 



20021 120061: iRarm'rez Velez fe Lopez Aristd [200 
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The results go in the direction of ubiquitous fields 
with strengths in the hG regime, which suppo r ts the 
Hanl e-effect conclusion of iTruiillo Bueno et al.l (|2004 
120011 . 

Recent results ([Martinez Gonzale z et al.ll2008a| ) indi- 
cate that there is not any apparent variation of the po- 
larimetric properties of the Fe 1 near-IR doublet at 1.56 
/Ltm with the heliocentric angle. This suggests that the 
magnetism of the quiet Sun remains approximately the 
same for every observing angle, thus giving weight to 
the fact that the distribution of magnetic field vectors is 
quasi-isotropic for the resolutions of ground based tele- 
scopes and discards a network-like scenario for the in- 
ternetwork magnetism. Based on these observations, 
lLopez Ariste et all (|2007f l suggest a scenario in which the 
observed Zeeman signals are a result of statistical fluc- 
tuations of a quasi-isotropic vector field. This scenario 
suggests a way to reconcile simultaneously the Zeeman 
and the Hanle results. It also gives an explanation to the 
reason why the unsigned magnetic flux seems to be the 
same regardless of the spatial resolution of the observa- 
tion. 

Recent observatio ns carried out with the spec- 
tropolarimeter (SP ; iLites et ail I2001D aboard Hinode 
(|Kosugi et al.l l2007t ) have led to the detection of large 
amounts of linear polarization signals in the Fe 1 lines 
at 630 nm at unprecedent spatial resolutions of ~0.3". 
The interpretation of these si g nals based on simple 
inversion methods (jLites et alJ l2008h or using more 
elaborate Milne-Eddingt on models (jOrozco Suarez et al.1 
l2007d; ILites et al.l 12008ft suggest the presence of mag- 
netic fields with a large horizontal component. The 
probability distribut i on of field inclinations given by 
lOrozco Suarez et"al"1 (|2007d ) suggests that the majority 
of fields found in the field-of-view are horizontal (with 
inclinations 7 = 90° ± 15°), with a certain (small) frac- 
tion of fields bein g almost vertical. It i s also n ecessary 
to point out that lOrozco Suarez et alJ (|2007al ) argued, 
with the aid of numerical magneto-hydrodynamical sim- 
ulations, that it is possible to get a good estimation of the 
magnetic field strength from Milne-Eddington inversions 
of Hinode data. 

If the apparent lack of variation of the po- 



to the average intensity) in the continuum surrounding the 630 nm 
total polarization normalized to the average intensity, J (Q 2 + U 2 + 
+ U 2 ) x / 2 d\/(I) > 0.0035 and J(V 2 )^ 2 dX/{I) > 0.003, respectively. 

larimetric properties w i th th e heliocentric angle 
((Martinez Gonzale z et al.l l2008af) and the apparent 
large presence of horizon t al fields at dis k center 
(lOrozco Suarez et alJ l2007d : ILites et alJ l2008| ) are put 
together, this would mean that, irrespective of the 
location in the solar disk, the fields would always be 
predominantly horizontal. This result indicates that 
the distribution of magnetic field vectors is indeed close 
to isotropic since an isotropic vector field distribution 
always presents a much larger amount of horizontal than 
vertical field (for all external observers). Therefore, 
it is of great interest to verify to what extent are the 
observations (with the presence of noise and the ensuing 
uncertainties) unambiguously pointing to predominantly 
horizontal fields. 

To this end, we apply Bayes-ME, the Bayesian in- 
ference code for Milne-Eddington a tmospheres devel- 
oped by lAsensio Ramos et alJ (l2007bh. to a sub-field of 
the lar ge fiel d-of-view presente d by lOrozco Suarez et all 
(|2007d ) and ILites et al.l (|2008h in order to analyze the 
effect of noise and degeneracies on the retrieved param- 
eters. 

2. OBSERVATIONS 

The observations are those analyzed by ILites et ail 
( 2008), which were obtained at disk center on 2007 March 
10 with the spectropolarimeter aboard Hinode. The an- 
alyzed map is a small sub-field of 7.4" x7.4" taken from 
a large field-of-view of 302 " x 162" . This small sub-field 
has been also analyzed by lOrozco Suarez et"ail |2007c) 
to investigate the properties of very quiet internetwork. 
The observed spectral region contains the Fe I doublet at 
630.1 and 630.2 nm with a resolution close to 3x 10 5 . The 
intensity image in the local continuum normalized to the 
average continuum intensity in the sub-field is shown on 
the left panel of Fig. [1] while the total polarization, de- 
fined as j(Q 2 + U 2 + V 2 ) 1/2 dX/{I), is shown in the right 
panel of the same figure. This figure also presents the 
regions of / (Q 2 + U 2 f/ 2 d\/{I) > 0.0035 in red contours 
and / (V 2 )^ 2 d\/(I) > 0.003 in blue contours. After cali- 
bration, the noise level amounts to 1.2x 10 _3 / c for Stokes 
V and to l.lxl0~ 3 / c for Stokes Q a nd U. For more de- 
tails on the observations, we refer to ILites et all (|2008l ). 

3. BAYESIAN INVERSION 
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3.1. Bayes-ME 

Although the Bayes-ME code has be en extensively 
described bv lAsensio" Ramos et all (|2007bD . we give here 
a brief overview because of its novel characteristics. The 
Bayes-ME code is an inversion code built under the 
fra mework of the Bayesian a pproach to inference (see 
e.g.. iNeall [19931 : lGregory||2005f ). Let M be a model that 
is proposed to explain an observed dataset D and let / be 
a set of sensible a-priori information about the problem 
(for instance, the spectral resolution of the observations). 
The model M is parameterized in terms of a vector of pa- 
rameters of length 7Vp ar , 9 G R^". Due to the presence 
of noise in the observations, any inversion procedure is 
not complete by just giving the values of the model pa- 
rameters that better fit the observations. The full solu- 
tion to the inference problem is to provide the posterior 
probability distribution function (pdf) p(9\D, I) that de- 
scribes the probability that a given set of parameters is 
compatible with the observables given the a-priori knowl- 
edge. As a consequence, statistically relevant informa- 
tion about one parameter (irrespective of the value taken 
by other parameters) can be obtained from this pdf by 
marginalization (integration) of the rest of parameters. 
In order to calculate the posterior distribution p(9\D, I), 
the Bayes theorem can be applied to obtain: 



p(0\D,I) 



p{D\e,i)p{e\i) 

p(D\I) ' 



(1) 



where p(D\I) is the so-called evidence, a quantity that 
in the context of parameter estimation is of no impor- 
tance (since it is a constant that does not depend on 
the model parameters 6) but that turns out to be cru- 
cial in the context of model selection. The distribution 
p{6\I) is denominated the prior distribution and contains 
all relevant a-priori information about the parameters of 
the model. Usually, unless some information is available 
about the value of some parameters, it is common to use 
uninformative (or vague) priors like bounded un iform dis- 
tributions or Jeffreys' priors (e.g., [Neaj|[l993). Finally, 
p(D\9, 1) is the likelihood, a distribution that character- 
izes how well a model with parameters 9 reproduces the 
observed dataset. The Bayes theorem states that the 
probability that a model M becomes plausible once the 
information encoded in the data D is taken into account 
depends on how plausible the model was without data 
and how well the model fits the data. 

Assuming that the observables are represented as a 
vector d £ R N , that the result of evaluating the model 
M at a given set of parameters 9 is the vector y € K w 
and that the observations are contaminated by a noise 
component characterized by the vector e e l w , we have: 

ck = yi + ei, Vi (2) 

When the chosen model parameters exactly correspond 
to those of the observed dataset, the distribution of dif- 
ferences yi — di has to follow the distribution of the noise. 
Assuming that the noise is Gaussian distributed with 
a standard deviation given by the vector ar € R , the 
likelihood function is, therefore, given by the following 
Gaussian distribution: 



p{D\9,I) oc ]Jexp 



(Vi -d l f 



2<if 



(3) 



In the case of Bayes-ME, the model M is based on 
the simplified Milne-Eddington (ME) description of stel- 
lar atmospheres. The emergent Stokes profiles are cal- 
culated as a linear combination of the emergent Stokes 
profiles obtained from an arbitrary number of compo- 
nents, each of them characterized by the standard ME 
parameters. Although the Bayesian approach will allow 
us in the future to automatically choose the optimal num- 
ber of components N c using model selection techniques, 
N c is fixed and chosen a-priori in the present version of 
BAYES-ME(in other words, iV c is part of the a-priori in- 
formation summarized by I). The vector 9 contains all 
the ME parameters over which wc do inference (the selec- 
tion of free parameters is completely configurable) , while 
the remaining parameters ME are set as fixed quanti- 
ties. The Bayes-ME code neglects the fixed parameters 
and they are not included in the vector 9, although the 
result is equivalent to the case in which the posterior 
distribution is marginalized (integrated) over them af- 
ter setting Dirac delta priors on them (see q3.2[) . The 
parameters of the ME atmosphere for each component 
are: the Doppler width of the line in wavelength units 
(AAdopp), the line-of-sight component of the macroscopic 
bulk velocity (i>los)> the gradient of the source function 
{fi), the ratio between the line and continuum absorp- 
tion coefficients (?7o), & line damping parameter (a) and 
the magnetic field vector parameterized by its strength, 
inclination and azimuth with respect to a given reference 
direction {B, 6b and \b, respectively). We refer, some- 
what arbitrarily, to AAdopp, vlos, ft, r]o an d a as the 
thermodynamical parameters (although their pure ther- 
modynamical character is quite diffuse). Additionally, 
we take into account the filling factor of each compo- 
nent, fi, subject to the constraint Yli=i fi = 1- The 
maximum number of free parameters is 9N C — 1. 

In order to sample the posterior distribution func- 
tion and carry out the ensuing marginalizations, Bayes- 
ME utilizes a Markov Chain M onte Carlo (MCMC; 
iMetropolis et al.lll953t lNea]||1993h scheme based on the 
Metropolis algorithm. The initial proposal density distri- 
bution is a multivariate Gaussian with diagonal covari- 
ance matrix that is set to 10% of the allowed range of 
variation of the parameters. After a configurable initial 
period, the proposal density is changed to a multivari- 
ate Gaussian with a covariance matrix that is estimated 
from the previous steps of the chain 1 . In order to im- 
prove convergence, the covariance matrix is multiplied 
by a quantity that assures that the acceptance rate of 
proposed models is close to 25%, a value t hat is the the- 
oretic al optimal value for simple problems (|Gelman et al.l 
[19911 . 

3.2. Marginalization 

In order to obtain the posterior probability distribution 
function for one parameter and give estimations and con- 
fidence intervals, we have to marginalize (integrate out) 
the rest of parameters: 

p(9i\D,I) = J d6 1 dO 2 ---de^ 1 dO l+1 ---d9 Np ^p{0\D,I). 

(4) 

1 This is an improvement over the version of Bayes-ME pre- 
sented by Ascnsio Ramos et al. (2007b) which used a proposal den- 
sity with diagonal covariance matrix. 
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(prior distribution): 



a,l= / P (e\D,i)iog 2 ^j§jp-de. 



(5) 



Fig. 2. — Histogram of Kullback-Leibler divergence between the 
posterior and prior distributions for the magnetic field strength 
(black line), inclination (red line) and azimuth (blue line). This 
indicates that the amount of information present in the observa- 
tions is sufficient to return posterior distributions clearly different 
from the prior distributions for the magnetic field strength and for 
the inclination. The information is heavily reduced for the azimuth 
of the field. 



Interestingly, the output of any MCMC code is a con- 
verged Markov chain for each parameter whose his- 
togram is just proportional to p(9i\D, I). Therefore, cal- 
culating marginal posterior distributions is just a matter 
of making histograms. Any relevant statistical informa- 
tion can be obtained from the distribution or summa- 
rized using, in the appropriate cases, the median and 
confidence intervals or upper/lower limits. 

3.3. This analysis 

Due to the general character of Bayes-ME, it is pos- 
sible to a dapt it to special inversion schemes like that 
applied by lOrozco Suarez et al.l (|2007cf ) for the inversion 
of Hinode quiet-Sun internetwork magnetic fields. We 
slightly modified Bayes-ME to suit the scheme followed 
by these authors since our aim is to analyze the relia- 
bility of their results. The model used consists of one 
magnetic Milne-Eddington component contaminated by 
the presence of a stray- light component. This stray- light 
contamination, characterized by a filling factor a is ob- 
tained locally as the average of the pixels within a box of 
1" centered in the pixel of interest. The magnetic filling 
factor (fraction of signal represented by a magnetic at- 
mosphere) is, therefore, 1 — a. In the analysis, we focus 
on the magnetic filling factor. 

A direct consequence of using a Bayesian inference ap- 
proach is that we do not need to establish a selection 
criterion to choose which pixels to invert. For those pix- 
els where the signal-to-noise level is too low, the ensuing 
posterior probability distributions will straightforwardly 
show the degeneracies present in the data and the re- 
sulting posterior distributions will automatically be very 
similar to the prior distributions. Therefore, we invert all 
the pixels in the field-of-view and select afterwards which 
pixels provide a sufficient amount of information to the 
inference process. Apart from the pure visual inspec- 
tion of the resulting marginal probability distributions, 
we use a criterium to characterize the pixels in which 
there is sufficient information. The criterium is based on 
the calculation of the Kul lback-Leibler divergence (KLD; 
iKullback &: Leiblerl I1951D between the two probability 
distributions p(8\D, I) (posterior distribution) and p(8\I) 



The KLD, always a positive quantity, measures the differ- 
ence in number of bits between transmitting samples of 
the distribution p(9 \ D, I) using a code based on p(9\D, I) 
and transmitting samples of the distribution p(9\D, I) us- 
ing a code based on the distribution p(9\I). The more 
different both distributions are, the larger the additional 
information required. Therefore, larger values of -Dkl 
indicate that both distributions are very different. We 
apply this measure to compare the final posterior prob- 
ability distribution after the data has been presented to 
the inference code and the prior distribution. When the 
data is informative, the posterior distribution is very dif- 
ferent from the prior distribution and -Dkl will be large. 
Since the priors for the stray-light factor, magnetic field 
strength, inclination and azimuth of the field are uniform 
distributions in the intervals [0, 1], [0, 1500] G, [0°, 180°] 
and [0°, 180°], respectively, the KLD simplifies to: 

/>oo 

D KL = log 2 (0 max -0 min )+ / p(8\D, I) log 2 p(6\D, I)d0, 



(6) 

where the last term represents the negative of the stan- 
dard definition of the entropy (sometimes called negen- 
tropy) of a continuous probability distribution function. 
This shows that the limit -Dkl = is reached when the 
posterior distribution is equal to the prior distribution 
and that it can potentially reach very high values for very 
informative posterior distributions (because the entropy 
is not bounded from above). According to the results 
presented below, pixels with Dkl ^ 0.5 can be consid- 
ered to have enough information to roughly constrain a 
given parameter, although a more conservative threshold 
would be Dkl > L 

4. RESULTS AND DISCUSSION 

After applying the Bayes-ME code to all the pix- 
els in the field-of-view, the results consist of converged 
Markov chains for each parameter. As an additional 
test for convergence, we rerun the code several times 
obtaining almost indistinguishable results. By making 
histograms of each chain, we obtain the posterior dis- 
tribution function marginalized over all parameters ex- 
cept the one of interest. These marginal posterior dis- 
tribution functions show how the information encoded 
on the observed Stokes parameters constrain every indi- 
vidual parameter. In principle, the correct answer to 
the inversion process would be to give such posterior 
distributions for each parameter and each pixel. Since 
this is obviously not feasible for presentation purposes, 
we follow a different approach and try to condense the 
important statistical information present in each chain. 
We focus only on the properti es of the magnetic field 
vector thanks to the results of [Wcstcndorp Plaz a et"aT1 
(1998), who demonstrated that the cross-talk between 
thcrmodynamical and magnetic parameters is of reduced 
importance in Milnc-Eddington inversions. The rea- 
son for this behavior has to be found on the large dif- 
ference between the response functions for the Milne- 
Eddington thermody namical and magnetic parameters, 
as elegantly shown bv lOrozco Suarez fc Del Toro Iniestal 
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(|2007f h Note that, due to the Bayesian approach that we 
follow, the marginal posterior distributions for the mag- 
netic field parameters automatically include information 
about all possible values of these thermodynamical pa- 
rameters weighted by their probability. 

4.1. Amount of information 

The first key point is to estimate how much informa- 
tion about every parameter can be inferred from the ob- 
servables. To this end, we analyze the Kullback-Leibler 
divergences. Figure [3] shows the histogram of Dkl for 
the marginal posterior distributions of the magnetic field 
strength, inclination and azimuth in all pixels. It is clear 
from this figure that the amount of information available 
to constrain the azimuth is very poor because Dkl is very 
close to zero for all the pixels. This is a consequence of 
the fact that almost all pixels present posterior distri- 
butions for the azimuth that are indistinguishable from 
the uniform prior distribution, as we demonstrate below. 
There is more information available for constraining the 
field strength and inclination (larger values of £>kl) al- 
though several pixels present posterior distributions that 
are still similar to the prior. Table [T] summarizes the 
percentage of pixels with relevant information for each 
parameter in the field-of-view. 

Of interest is also the analysis of the horizontal vari- 
ation of £>kl shown in Fig. [3] The results clearly indi- 
cate what has been summarized in Fig. O For clarity, 
we have also indicated with contours the regions where 
Ic/{Ic) = 1-02, which may serve as a rough indication 
of the position of the brighter parts of granules. Several 
interesting points can be inferred from these plots. 

First, a large amount of pixels present large values of 
-Dkl(1 — ot) (bright points in the upper left panel of the 
figure). These pixels have, as we show later in Fig. QJ 
posterior distributions for the magnetic filling factor with 
a very conspicuous peak close to zero. They correspond 
to pixels in which the polarization signal is extremely low 
and well below the noise. In such a situation, and since 
the intensity profile is non-zero, the best fit is accom- 
plished with ^100% stray-light contamination. For this 
reason, the very same white pixels for the -Dkl(1 — oc) 
map are typically dark in the -Dkl(-B) map because the 
noise rapidly destroys the information needed to give a 
reliable value of the magnetic field strength. It is impor- 
tant to point out that these results are only valid for the 
present noise level of the observations and that results 
would surely change if a better signal-to-noise ratio had 
been achieved. 

Second, it is obvious from the middle panels of Fig. 
[3] that there is less information to constrain the inclina- 
tion and azimuth of the field (especially for the azimuth) . 
However, it is interesting to point out that there is appar- 
ently more information available in intergranular lanes 
than in granules, in many cases corresponding to pixels 
with the largest polarization signals. These polarization 
signals in the integranular lanes are mainly circular po- 
larization, as shown in Fig. [TJ It is also obvious from Fig. 
[3] the appearance of patches in the field-of-view where 
the field strength, inclination and azimuth can be nicely 
constrained. Note that, according to Fig. [TJ they corre- 
spond to patches where the linear polarization signal is 
large. This demonstrates, as we shown in f|5]and Fig. [5J 
the well-known fact that a linear polarization signal well 



TABLE 1 

Pixels with Kullback-Leibler divergences above different 
thresholds. 



Parameter 


N(D KL ) > 0.5 


N(D KL ) > 1 


N(D KL ) > 2.5 


a 


100.0% 


100.0% 


77.5% 


13 


64.8% 


36.9% 


9.2 % 


Ob 


49.8% 


33.7% 


1.5% 


Xb 


7.0% 


3.8% 


0.8% 


B \\ 


90.6% 


61.0% 


14.1% 





100.0% 


100.0% 


99.7% 



above the noise level produces a well constrained inferred 
magnetic field vector. Particularly frustrating is the fact 
that, at this noise level and apart from these patches, the 
information for constraining the azimuth of the field is 
almost negligible. 

Finally, the lower panels present the horizontal varia- 
tion of the Kullback-Leibler divergence for B» , the longi- 
tudinal component of the magnetic field vector (lower left 
panel) and (f> = (1 —a)Bu, the longitudinal magnetic flux 
density (lower right panel) . The information available for 
the linc-of-sight component of the magnetic field vector 
is very similar to that of the field itself but slightly larger 
because it is augmented with the knowledge of the incli- 
nation angle. On the contrary, the longitudinal magnetic 
flux density shows a very clear correlation with the gran- 
ulation pattern. All values of -Dkl(</>) are clearly larger 
than 1, showing that there is sufficient information in 
all the pixels to obtain information about this quantity, 
which clearly changes smoothly from pixel to pixel. 

4.2. Marginal posteriors 

In order to understand the typical values of -Dkl for 
which we can consider that the posterior distribution is 
clearly different from the prior distribution, we show in 
Fig. [4] examples of the posterior distributions for the 
magnetic filling factor, magnetic field strength, inclina- 
tion and azimuth for different values of £>kl- We also 
present the posterior distributions for the longitudinal 
component of the field and for the longitudinal mag- 
netic flux density. They show, in general, that nicely 
constrained parameters present £>kl 1- in order to 
facilitate the comparison, we also overplot the prior dis- 
tribution in dashed lines. 

As already discussed, the stray-light contamination is 
usually well constrained, at least for values of Z?kl(1 — 
a) > 1.3 (see upper left panel of Fig. |4]). Smaller val- 
ues present posterior distributions that can be considered 
commensurate with the prior distribution. Intermediate 
values typically correspond to those points in which the 
stray-light contamination amounts to a ~ 0.8 on aver- 
age, while the Bayesian results for those points with large 
values of -Dkl(1— a) point to an al most non-magnetic at- 
mosphe re. As already discussed bv lOrozco Suarez et al.l 
(2007b]), this behavior is a consequence of the fact that 
the stray-light contamination is mainly sensitive to the 
Stokes I profile and can be constrained fairly well irre- 
spective of the noise level in the rest of Stokes parame- 
ters. 

A similar behavior is seen in the marginal posteriors for 
the magnetic field strength (upper right panel of Fig. [3}, 
although the results are clearly less informative. When 
Dkl(B) < 0.1, we can state that the posterior distri- 
bution is very similar to the prior distribution (marked 
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Fig. 3. — Horizontal variation of the Kullback-Leibler divergence between the posterior and prior distributions for the magnetic filling 
factor (upper left panel), the magnetic field strength (upper right panel), field inclination (middle left panel), field azimuth (middle right 
panel), longitudinal component of the field (lower left panel) and longitudinal magnetic flux density (lower right panel). The contours 
indicate regions where I C /(I C ) = 1.02. 
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Fig. 4. — Marginal posterior distribution functions for the magnetic filling factor (upper left panel), magnetic field strength (upper right 
panel), field inclination (middle left panel), field azimuth (middle right panel), longitudinal component of the field vector (lower left panel) 
and longitudinal magnetic flux density(lower right panel) for several pixels with different values of the Kullback-Leibler divergence. The 
selected pixels are representative of zones where there is no information about the parameters (small values of Dkl) and those in which 
the information is enough to restrict to some extent their values (large values of -Dkl)- The posterior distributions have been smoothed 
with a Gaussian kernel for aesthetic purposes. For comparison, the prior distribution is shown in dashed lines. 



as a horizontal dashed line), so that these pixels do not 
contain enough information for constraining this param- 
eter. When Dkl(B) increases, the posterior distribution 
functions usually show the same behavior, peaking at 
magnetic fields close to zero but with an extended tail 
towards higher values. For such distributions with ex- 
tended tails, the only quantity with statistical meaning 
that we can give to summarize the results is an upper 
limit (with 68% or 95% confidence). Although an anal- 
ysis of these upper limits will be presented in M4. 31 it 
is possible to sec from the upper right panel of Fig. [4] 
that they are systematically favouring hG fields with re- 
spect to stronger fields. The set of points associated to 
the bright patch of the middle panels of Fig. [3] pro- 
duce marginal posterior distributions with a Gaussian- 
like shape peaking in the 400-600 G range but with tails 
that completely discard fields in the kG regime. This 
represents one of the key conclu sions of our work. 

After the claims raised by iLites et ail (|2008f) and 
lOrozco Suarez et al.l (|2007cf ) about the large amount of 
very inclined magnetic fields (with a distribution that 
peaks near 90°), it is interesting to analyze how much 
information is still encoded in the Stokes profiles of the 
Hinode observations under study as analyzed with the 
suggested Milne-Eddington model. It is important to 
point out that, although the pres ence of a la rge amount 



of inclined fields was inferre d by ILites et alj (|2008{ ) and 
lOrozco Suarez et"al"1 (|2007d ) with data with the same 
no ise level as those p resented here, it was also confirmed 
bv lLites et al.l (|2008 ) with observations with a noise level 
a factor ~ 4 smaller. The middle left panel of Fig. 
[4] shows the marginal posterior probability distribution 
functions for the field inclination. We have verified that 
it is fundamental to have -Dkl($b) > 0.5 in order to end 
up with a posterior distribution function clearly different 
from the prior distribution (dashed line). According to 
Table [T] and Fig. [2 the number of points fulfilling this 
condition is close to 50%. We have to point out that 
many of these points present posterior marginal distri- 
butions similar to the green or blue curves in the middle 
left panel of Fig. 2] Such posteriors, that discard half of 
the range of inclinations but give equal probability to the 
other half, are a consequence of the small Stokes V signal 
that is still enough to univoquely select the direction of 
the field but not its inclination. In other words, the po- 
larity of the field is easily recovered but it remains impos- 
sible to constrain more the inclination. This means that, 
even for pixels with low signals (well below the thresh- 
old), there is enough information to, at least, detect the 
polarity of the field. Obviously, when the filling fraction 
of the magnetic component is very small and the mag- 
netic field is in the weak field regime, the field cannot be 
recovered uniquely. This translates into posterior distri- 
bution functions with long tails for the field strength and 
the field inclination, as a conseque nce of the marginaliza- 
tion ov er the rest of variables (see lAsensio Ramos et al.1 
I2007bh . 

Concerning the field azimuth (middle right panel of 
Fig. 3]) , the results are not very encouraging. The ma- 
jority of marginal posterior distributions commensurate 
with the prior distribution. Only in approximately 1% 
of the field-of-view we find marginal posteriors with a 
clear peak favoring a certain subset of azimuths (those 
corresponding to the white patches of the middle right 



panel of Fig. [3]), while in ~4% of the points, the infor- 
mation encoded in the Stokes profiles allows to roughly 
determine the azimuth but with a very large error bar. 

The lower two panels of Fig. [4]present several posterior 
distributions for Bu and for pixels with different val- 
ues of the Kullback-Leibler divergence. The dashed lines 
show the prior distributions, whose functional forms are, 
assuming uniform priors for a, B and 9b- 
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for the longitudinal magnetic flux density. Note that 
p(Bu = 0) and p(4> = 0) tend to infinity in this case 
although the distribution has finite area. We believe it 
is important to point out that the product of random 
variables distributed uniformly can result in highly non- 
trivial distributions. Consequently, it is extremely im- 
portant to understand which is the inherent prior infor- 
mation introduced in any inversion method because, when 
the data is non-informative, the parameters recovered can 
be strongly influenced by the prior. 

The posterior distributions for B» for small values of 
Dkl(B\\) resemble the prior distribution. They become 
much more distinct when enough information is encoded 
in the Stokes profiles, thus resulting in peaked posterior 
distributions without extended tails. Interestingly, as a 
consequence of the multiplication by the well-determined 
magnetic filling factor, all the posterior distributions for 
the longitudinal flux density are quite dissimilar to the 
prior, except perhaps those points with small values of 
the Kullback-Leibler divergence, i.e., -Dkl(</>) ~ 1-7. 

4.3. Parameter estimation 

According to the results of the previous section, the 
magnetic filling factor and the upper limit of the mag- 
netic field strength can be nicely inferred for a relatively 
large amount of pixels. Likewise, it is possible to set 
constraints on the field inclination for a subset of the 
field-of-view, while almost no pixels contain information 
about the field azimuth. Since the posterior distributions 
for I — a are fairly peaked, we summarize them using the 
median value, together with the upper 95% confidence 
limits. It is calculated as the value that encloses 95% of 
the total area of the posterior distribution starting from 
1 — a = 0. The results are shown in the upper left panel 
of Fig. where we plot the histogram of median values 
in black and the histogram of the 2a limit in red, respec- 
tively. The vertical dashed line ind icates the position of 
the mo st abundant value found by lOrozco Suarez et al.l 
(|2007d) . Note that if the plot is done in linear scale 
almost all the points collapse to I — a « 0. Concern- 
ing the confidence intervals, the histograms indicate that 
they are relatively narrow, completely discarding mag- 
netic filling factors with I — a > 0.4 at 95% confidence. 
This upper l imit is in rough agreement w ith the results 
presented bv lOrozco Suarez et al.l (|2007d ). 

Because the shape of the posterior distributions for the 
magnetic field strength only allows us to put upper limits 
in this parameter, we show in the upper right panel of 
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Fig. 5. — Histograms showing the parameters estimated from the marginal posterior distribution functions for the magnetic filling factor 
(upper left panel), upper limit of the magnetic field strength (upper right panel), inclination (lower left panel) and longitudinal magnetic 
flux density (lower right panel) for those pixels containing enough information. Note that this include 100% of the pixels for I — a and 
4>, 36.9% of the pixels for B and only 33.7% of the points for 8b- For the magnetic field strength, the results have been separated into 
granules and intergranules. 

Fig. [5] the histograms of the horizontal variation of the 
68% and 95% upper limits of the distributions only for 
those pixels having Dkl(B) > 1 (roughly one third of the 
field-of-view, according to the results shown in Table [T]). 
For illustration, we have also divided the field-of-view 
into pixels belonging to granules and to intergranular 
lanes. The selection criterion we have used is very sim- 
plistic and it is based on selecting pixels with I C /(I C ) > 
1.02 as granules and pixels with l c j (I c ) < 0.98 as inter- 
granules. The histograms indicate a quite unclear dif- 
ference between the limit field in "granules" and "inter- 
granules" , with the histogram of the 68% confidence limit 
field around 200 G for both cases. The histograms are ob- 
viously shifted towards higher field strengths for the 95% 
confidence histograms but the results undeniably discard 
fields above ^900 G. According to the information avail- 
able in the Stokes profiles observed with Hinode, these 
histograms are the only statistically relevant result that 
one can give about the magnetic field strength, except 
for some individual points where the posterior distribu- 
tion presents a clear Gaussian-like shape. It is important 
to note that these results do not discard probability dis- 
tribution functions peaking at muc h smaller fields lik e 
those present in MHD simulations (jVogler et alJ 1 20051 . 
or those inferred from infrared data (e.g., iLinl 119951 : 



Khomenko et all 120031: iMartmez Gonzalez et~aH l2006ri 
Martinez Gonza lez et al. 2008 b) or lines with hyper- 
fine structure ([Ramirez Velez et alJ 120081 ). However, 
they discard pixels in the field-of-view having fields 
above 900 G. In this sens e, the results obtained by 
lOrozco Suarez et"aTI (|2007ch using standard Levenberg- 
Marquardt inversion schemes are inside the limits posed 
by our Bayesian analysis. 

Only around 50% of the Stokes profiles observed in the 
field-of-view contain enough information for introducing 
constraints on the inclination (£>kl(#b) > 0.5). If we 
summarize the inferred inclination by the median value 
of the posterior distribution (without taking into account 
its spread), the result is shown in the lower left panel of 
Fig. [3 Because of the typically extended shape of the 
posterior distributions for the field inclination, the results 
are therefore of reduced interest as representative of the 
statistical properties of the magnetic field inclination in 
the whole field-of-view. The histograms of 68% confi- 
dence intervals are also presented. The results indicate 
that the median values peak close to 30° and 140° with 
a dispersion around these values. This is a consequence 
of the fact that for many points we can only distinguish 
the polarity of the field (with the posterior distribution 
resembling, for instance, the curve for Dkl{Qb) = 1-1 in 
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Fig. 6. — This figure shows the histogram of maximum ampli- 
tudes of Stokes Q, U and V. The histogram in dashed line is the 
corresponding to pure noise calculated using order statistics. The 
vertical dashed line indicates the threshold of 4.5 times the noise 
level. 

the middle left panel of Fig. and the median value is 
close to the center of the [0°, 90°] or [90°, 180°] intervals. 

Finally, the lower right panel of Fig. [5] represents the 
median (in black line) and the ±ler values (red and blue 
lines) of the magnetic flux density as obtained from the 
marginal posterior distributions of the lower right panel 
of Fig. |H For comparison, we also show the histogram 
of the magnetic flux density obtained under the assump- 
tion of the weak field regime for those pixels with signals 
above 4.5 times the noise level. The magnetic flux den- 
sity histogram estimated under the assumption of weak 
field regime lies inside the Bayesian ±lcr confidence re- 
gion. This is a consequence of the fact that, even if the 
cross-talk between the strength and inclination of the 
magnetic field vector and the magnetic filling factor can- 
not be disentangled, the product <fi = (l — a)B cosOb can 
be estimated with confidence. 

4.4. Available information 

The fundamental point to clarify at the light of the re- 
sults shown in the previous sections is which is the mech- 
anism (or combination of mechanisms) that is destroying 
the information about the magnetic field vector in the ob- 
served Stokes profiles. To this end, we first analyze the 
fraction of pixels with signal above the noise level. Fig- 
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Fig. 7. — Inferred magnetic filling factor, the upper limit of the 
field and the inclination of the field versus the maximum amplitude 
of the polarization profiles, i.e., max{Q(A)// c , U(X)/I C , V(X)/I C }. 
The dashed horizontal line shows the noise level and the solid line 
presents the threshold of 4.5 times the noise level. This shows that 
pixels with amplitudes below and above the threshold of 4.5 times 
the noise level present similar behaviors. 



ure[6]shows the histogram of the maximum amplitude of 
Stokes Q, U and V in each pixel (the quantity chosen by 
lOrozco Suarez et"aT1l2007ci to set their threshold). The 
histogram in dashed line presents what one would obtain 
for a pure noisy case with the same standard deviation 
of the Hinode data, while the vertical dashed line indi- 
cates the threshold of 4.5 times the noise level. Note 
that the probability distribution function of a random 
variable that results from taking the maximum of a fi- 
nite number of normal random quantities (the pixels in 
the wavelength direction of the ca mera) can be obtained 
using techniques of order statistics fbavid 198lJ). For the 
case of Gaussian noise and 90 points in wavelength direc- 
tion, this gives the skewed distribution shown in Fig. [Sj 
We have empirically verified that the distribution peaks 
at ~ 2.5c, where a is the standard deviation of the noise. 

In it clear that Stokes Q and U profiles arc mainly 
dominated by noise in a large fraction of the pixels, with 
a tail towards large amplitudes that is slightly heavier 
than for the case of pure Gaussian noise . The thresh- 
old chosen bv lOrozco Suarez et al . (2007c) appropriately 
isolates these points. On the contrary, the Stokes V am- 
plitudes are not dominated by noise for amplitudes above 
the threshold. This means that, in the majority of the 
pixels, only the (possibly modified) information encoded 
in Stokes / and Stokes V remains. Note that the infer- 
ence is carried out using the full wavelength variation of 
the profile, so that it is certain that, even if the ampli- 
tude is close to the noise level, some information may still 
remain in the correlation between different wavelength 
points along the profile. When the stray-light contami- 
nation is large in the pixels, the information present in 
Stokes V is also efficiently masked by the noise because 
the circular polarization signal decreases. Consequently, 
once the stray-light contamination is set by adapting the 
model until fitting Stokes /, the upper limit in the field 
strength is essentially set by the remaining amplitude of 
Stokes V. 

The fact that noise is not the only contributor to 
the lack of information in the Stokes profiles can be 
reinforced thanks to the results presented in Fig. [7) 
In this figure we plot the median values of the mag- 
netic filling factor, the 68% upper limit of the magnetic 
field strength and the median value of the inclination 
of the field for those pixels in which £>kl(1 — o) > 1, 
Dkl(B) > 1 and Dkl{@b) > 0.5 versus the value 
of max{Q(X)/I c ,U(X)/I c ,V{X)/I c }. We also show the 
noise level as a horizontal dashed line and 4.5 times the 
noise level as a horizontal solid line. All pixels present 
significant signal of Stokes Q, U or V, with this signal 
always above ~ 2a. Except for the stray-light contami- 
nation, the points above and below the threshold set by 
lOrozco Suarez et"al"1 (|2007ch present roughly similar re- 
sults and seem to be independent of the amplitude of the 
profiles. A different behavior is found for the stray-light 
contamination. Profiles with signals below the thresh- 
old show systematically very large stray-light contamina- 
tions (with 1— a ~ 0.01), while pixels above the threshold 
show systematically smaller contaminations. 

Although apparently obvious, another quantitative 
conclusion that we can extract from our analysis is that 
one needs to detect signal in Stokes Q and U in order to 
get a good estimation of the inclination. This is demon- 
strated in the right panel of Fig. [8] We present the 
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Fig. 8. — Posterior distributions for the field inclination in two 
pixels with sizable Stokes Q and U signals (black solid and dashed 
lines) and in two pixels with negligible Stokes Q and U signals (red 
solid and dashed lines). 

posterior distribution function for two pixels with neg- 
ligible Stokes Q and U that are not detected above the 
noise level (red curves). The posterior distributions agree 
with the prior uniform distribution, so no information 
about the inclination is available. On the contrary, for 
the two pixels for which the Stokes Q and U signal are 
well above the noise level the posterior distribution func- 
tions are clearly different from the prior distribution. 

5. QUASI-ISOTROPY OF THE MAGNETIC FIELD 

An obvious question that arises after our analysis is 
whether it is possible to extract some solid conclusions 
about the field inclination from the information available 
in the observations. In order to answer to this question, 
we have calculated, for each pixel, the following quanti- 
ties, that give the fraction of the total area of the pos- 
terior distribution function that is enclosed in a given 
interval: 

f(6 1 ,6 2 )= p{e B )&9 B , g(xi,X2)= p(XB)d XB . 
Je 1 J X i 

(9) 

The posterior marginal distributions p(6b) and p(xb) 
are normalized to unit area. The quantity / represents 
the probability that the inclination of the field lies be- 
tween 9i and 02- In an isotropic field, one finds the same 
amount of vectors whose inclination lies between 0° and 
60° (or between 120° and 180°) than between 60° and 
90° (or between 90° and 120°). In other words, using Eq. 
®, /(|/i| < 0.5) = f(\n\ > 0.5) = 1/2, with n = cos6> B . 
Correspondingly, the quantity g represents the proba- 
bility that the azimuth of the field lies between xi an d 
X2- In an isotropic field, one finds the same amount of 
vectors with azimuths between 0° and 45° than between 
45° and 90° (the same happens for azimuths between 
90° and 180°). According to our imposed prior on the 
azimuth, we are not able to distinguish a purely ran- 
dom azimuth from a situation in which we have a non- 
constrained azimuth. The left panel of figure [9] shows the 
histogram of the quantity f(\n\ < 0.5) for those points 
whose polarimetric signal is above the threshold (dashed 
line) and below the threshold (solid line). Note that, in 
an isotropic field, the histogram of f(\p,\ < 0.5) should 
peak at 1/2 (with some natural dispersion). Although 
the histogram of pixels with signals below the thresh- 
old (dashed line) clearly point towards a quasi-isotropic 



distribution, the histogram of pixels with signals above 
the threshold (solid line) indicates a clearly non-isotropic 
distribution. In this last case, the histogram of high- 
inclination fields (fields with < 0.5) peak at values 
conspicuously smaller than 1/2. For completitude, the 
histogram for low-inclination fields (fields with \fi\ > 0.5) 
peaks therefore at values larger than 1/2. Note that, 
according to our imposed prior (uniform in 9b), in the 
complete absence of information on the Stokes profiles, 
/(|/Lt| < 0.5) should converge to the value 1/3. The fact 
that the histogram peaks close to 1/2 means that infor- 
mation about the quasi-isotropicity of the magnetic field 
vector is still available. Concerning the azimuth, we plot 
in the right panel of Fig. [5] the histogram of the pixel- 
to-pixel variation of the posterior mass for the azimuth 
between 45° and 135°. According to the previous dis- 
cussion, this should be close to 1/2 for an isotropic field, 
although this should also be the case for completely non- 
constrained azimuths. We find that this is the case in 
general for all points in the field-of-view, irrespective of 
the amplitude of the polarimetric signal, thus suggest- 
ing the idea that the azimuth of the field is not clearly 
organized. Therefore, for pixels with polarimetric signal 
below 4.5 times the noise level, the field is apparently 
quasi-isotropic (as a vector field). We point out that this 
result is in agreement with the quasi-isotropic characte r 
of the field found by iMartmez Gonzalez et al.l ((2008a) . 
A more deep analysis of improved data is mandatory to 
evaluate this property of the field. 

6. COMPARISON WITH x 2 -MINIMIZATION 

The model we propose for explaining the obs e rvation s 
is exactly that used by lOrozco Suarez etlrtl (|2007d) . 
Since we are using uniform priors for all the variables 
and we are sampling the full posterior distribution, we 
are also capable of recovering their results by just looking 
at the combination of parameters that gives the maxi- 
mum value of p(9\D,I). Such model is also the global 
minimum of the X m erit function and should be e quiva- 
lent to that found bv lOrozco Suarez et alj (|2007cD using 
Levenberg-Marquardt techniques if the method was suc- 
cessful in locating the global minimum. We point out 
that the maximum-a-posterior (MAP) model, although 
it represents the model that better fits the observations, 
is statistically equivalent to all those models inside the 
68% confidence region, as stated above. Consequently, 
the presence of noise induces that the MAP model can- 
not be preferred over the rest of compatible models based 
on statistical justifications. 

The correct Bayesian way of giving constraints for 
individual parameters is done after marginalizing the 
rest of parameters (which are accounted for with their 
associated probabilities). Unless the multidimensional 
posterior distribution is peaked (Gaussian-like, for in- 
stance), the MAP model and the marginal value for 
each parameter do not coincide. This is especially rel- 
evant when parameters are highly degenerate, which is 
the case with this dataset. We have verified that the 
M AP results agree well with t hose presented in Fig. 2 
of lOrozco Suarez et al.l (|2007cf ). as shown in Fig. [TU1 
but they usually differ from the value at the peak of the 
marginal distribution. This is a clear indication of the 
presence of strong degeneracies. All points whose polari- 
metric signal is not above 4.5 times the noise level are set 
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Fig. 9. — Histogram of the posterior distribution mass for the inclination normalized to the total area of the posterior (left panel) in 
different inclination regimes. The right panel shows the posterior distribution mass for the azimuth between 45° and 135°. The results 
suggest that the field is quasi-isotropic for those points whose polarimetric signal is below 4.5 times the noise level. 



to black in Fig. [TO] The histogram of MAP field inclina- 
tions presents a shape very similar to that found in the 
lower left panel of Fig. [5j with two peaks symmetrically 
placed around 90° . Interestingly, the MAP estimation of 
the field inclination for those pixels with amplitude sig- 
nals below the threshold tends to be very close to 90°. 
Indeed, these pixels also correspond to the points with 
Dkl{0b) < 0.5. T his suggests that many of th e points 
giving 9 B ~ 90° in lOrozco Suarez et all (|2007d) should 
be reconsidered within a Bayesian framework in order to 
investigate whether this value is really representative or 
they are just a subproduct of the lack of information on 
the Stokes profiles. The fact that many of these points 
with inclination angles of the magnetic field close to 90° 
are located inside regions with low polarization ampli- 
tud es and/ or surrounding region s of large signal (see Fig. 
1 of lOrozco Suarez et al.ll2007cT ) points towards the sec- 
ond possibility. This analysis is left for the future. 

7. CONCLUSIONS 

We have applied the Bayes-ME inference code to 
Hinode spectropolarimetric observations of the internet- 
work quiet Sun. The deep analysis carried out under 
the Bayesian inference approach and using a simplified 
Milne-Eddington model of one magnetic component with 
a stray-light contamination leads us to conclude that 
there is information in the observables to put upper lim- 
its to the magnetic field strength and to give a good 
estimation of the stray-light contamination. The results 
indicate that the magnetic field strength is clearly in the 
hG regime with 95% confidence, so that fields in the kG 
regime are effectively discarded. All the results presented 
in this paper depend on the assumption that the model 
proposed explains the observables. It might be the case 
that a simpler or more complex model allows us to ex- 
tract more information from the observables because it 
inherently reduces some of the degeneracies. This is an 
issue left for future investigation. 



Concerning the field inclination, the i nformation 
presen t in t he Hinode Stokes profiles used bv lLites et al.l 
(|2008ft and lOrozco Suarez et all (|2007d ) is clearly not 
enough to tightly constrain it. However, since after the 
Bayesian analysis we have in hands the full posterior 
distribution for each parameter and each pixel, we have 
been able to conclude that the amount of high- inclination 
fields (with 9b < 60° and 0b > 120°) and low- inclination 
fields (with 60° < 9b < 120°) is essentially the same for 
pixels with polarimetric amplitudes below 4.5 times the 
noise level. This suggests that the magnetic fields of part 
of the quiet solar photosphere are well described by an 
quasi-isotropic distribution. This supports the results of 
iMartmez Gonzalez et al . (2008a). For those pixels with 
polarimetric amplitudes above the threshold, the incli- 
nations are clearly constrained. 

This paper also demonstrates that, although Bayes- 
ME 2 is computationally more intensive than standard 
Levenberg-Marquardt codes, it is feasible to analyze real 
data under the Bayesian framework. The analysis of the 
~ 2500 pixels of the map was done in less than 8 hours in 
a standard desktop computer, at a rate of ~10 seconds 
per pixel. In spite of the increase in computing time, the 
amount of information inferred from the analysis is far 
more complete than that from standard gradient descent 
inversions. This has allowed us to verify which parame- 
ters are constrained by the observations and permits the 
user to yield conclusions based purely on the observa- 
tions, clearly showing when conclusions can be flawed by 
a-priori assumptions. 
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